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i-H ' Context. Statistical properties of H n region populations in disk galaxies yield important clues to the physics of massive star formation. 

Aims. We present a set of H n region catalogues and luminosity functions for a sample of 56 spiral galaxies in order to derive the most 
general form of their luminosity function. 
t-H ■ Methods. H n region luminosity functions are derived for individual galaxies which, after photometric calibration, are summed to 

' form a total luminosity function comprising 17,797 Hn regions from 53 galaxies. 

£ — ' Results. The total luminosity function, above its lower limit of completeness, is clearly best fitted by a double power law with a 

significantly steeper slope for the high luminosity portion of the function. This change of slope has been reported in the literature for 
\Q ■ individual galaxies, and occurs at a luminosity of log L = 38.6 + 0. 1 (L in erg s~' ) which has been termed the Stromgren luminosity. A 

^\ • steep fall off in the luminosity function above log L = 40 is also noted, and is related to an upper limit to the luminosities of underlying 

massive stellar clusters. Detailed data are presented for the individual sample galaxies. 
\Q , Conclusions. The luminosity functions of H n regions in spiral galaxies show a two slope power law behaviour, with a significantly 

steeper slope for the high luminosity branch. This can be modelled by assuming that the high luminosity regions are density bounded, 
though the scenario is complicated by the inhomogeneity of the ionized interstellar medium. The break, irrespective of its origin, is 
of potential use as a distance indicator for disc galaxies. 
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C\j 1. Introduction bers of high luminosity Hn regions to test the suggestion that 

the break luminosity shows low scatter. In three articles (Rozas, 
In their paper giving a major overview of the luminosity func- Beckman & Knapen 1996; Rozas et al. 1999; and Rozas, Zurita 



tions (LFs) in Ha emission of the H n regions in disk galaxies & Beckman 2000), eight galaxies were measured and all showed 

^ . Kennicutt, Edgar & Hodge (1989; hereafter KEH) noted that type n LFs , most strikingly NGC 7479, which has the largest 

03 an important subset of their objects showed LFs with a clear number of luminous H n regions 

break in slope at a luminosity of L(Ha) = 38.7 dex (units of In spite of the general result reported in KEH, and in spite 
ergs 1 ). At lower luminosities the power law slope is flatter, of the type n LFs found in individual objects, the question has 
whereas above the break it is steeper. KEH called those LFs been raised persistently whether the break is a real feature or an 
showing a break "type II", suggesting that any galaxy will ex- Two types of doubts can be found in the literature. One 
hibit type II behaviour if it has sufficient high luminosity H n claims that though breaks can be found, they are not a universal 
regions. A galaxy with sufficient luminous H n regions to show physical f ea t U re, as the break luminosity varies over quite a wide 
a clear break, M51, was observed by Rand (1992), who mea- range between objects. One example here is by Thilker, Braun 
sured it at L(Ua) = 38.6 dex. Rand's work, using CCD data, was & Walterbos (2000) who derived an LF for M5 1 using an auto- 
especially precise, and revealed, as well as the break, a narrow mated met hod, finding a break at 38.9 dex, compared with the 
peak in the LF close to the break luminosity. In the LF of M33, value of 38 . 6 by Rand ( i 992 ) who used an interactive, region by 
Hodge et al. ( 1 999) found no break, which is not surprising since region met hod. The point raised is whether for statistical reasons 
the function is defined only up to L(Ha) = 38.4 dex, with only a conne cted with the chosen binning parameters a clear break lu- 
dozen regions above L(Ha) = 38 dex. Further examples of this minosity cannot be derived, given the low number of regions per 
type are NGC 6822 and the dwarf galaxy Holmberg II (Hodge, luminosity bin. Others have not found consistent evidence for 
Strobel & Kennicutt 1994), smaller galaxies with small Hn re- breaks, and certainly not at a specific luminosity, among these 
gion numbers at L(Ha) > 38 dex. On the other hand, M101 has b eing Gonzalez Delgado & Perez (1997) in an Ha survey of 27 
some 60 Hn regions with L(Ha) > 38.5 dex and shows a clear ga i axies selected for their nuclear activity. Their sample corn- 
break at L(Ua) 38.6 dex (Scowen, Dufour & Hester 1992). Our prised SO me 2000 Hn regions, i.e., * 75 per galaxy in a lumi- 
group has published LFs for galaxies selected to have large num- nosity range from L(Ha) = 37 dex to 39.5 dex. The numbers per 
bin above 38.5 dex are on average quite small for each galaxy 
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On completing a recent imaging survey in Ha of a set of 
nearby galaxies (Knapen et al. 2004) we saw that we could make 
a significant test of the suggestion that a dual slope LF with a 
well defined break gives a better general description than a single 
power law slope, and that the original result of KEH, with its 
clean break at a specific luminosity, is present in a data set with 
a sufficient statistical base. Our observations have allowed us to 
isolate and catalogue almost 18,000 Hn regions in 56 galaxies. 
Here we give the key results of this study and also offer on-line 
catalogues for the full set of objects. 

2. Sample selection and Ho imaging 

2.1. Sample selection and observations 

We have used the data set published by Knapen et al. (2004), 
which contains a full set of continuum-subtracted Ha images of 
57 relatively face-on (i <50°), nearby, Northern, spiral galaxies. 
Most of the images were obtained with the 1 m Jacobus Kapteyn 
Telescope (JKT), but some with other telescopes or from the lit- 
erature. Most of the Ha images were taken through « 50A wide 
filters matched to the recessional velocity of each galaxy. Full 
details on the observations are given in Knapen et al. (2004), 
whereas the sample galaxies are listed, and their distances, plate 
scales, and resolutions given, in Table 1 (online only). 

2.2. The effects of seeing on Hn region selection 

The FWHM seeing in our Ha images varies from 0.8 to 
3.7arcsec (1.53 arcsec in the median), which corresponds to a 
range in spatial resolution of 30 to 295 pc (137pc in the me- 
dian). We thus need to consider possible effects of blending 
on the resulting Hn region catalogues. The two major effects 
are that fainter Hn regions may be spatially coincident with 
larger, brighter, H n regions, reducing the number of faint H n 
regions observed, and overlapping of Hn regions, which may 
lead to a measured luminosity higher than the true luminosity. 
KEH showed that whereas degrading images up to resolutions 
of 200 pc had little effect on the luminosities of interest here, 
and that degrading them to a resolution of 200-400 pc will af- 
fect the faint end of the luminosity function while preserving the 
shape of the upper LF, blending seriously affects the entire LF 
and causes increases in the luminosities of first ranked Hn re- 
gions in images with a resolution of 300-500 pc. Rand (1992) 
modelled the main effects of blending for the arm LF of M51, 
and found that they did not significantly affect the shape of the 
LF. 

More recently, Ha LFs have been derived from HST imaging 
(Pleuss, Heller & Fricke 2000; Scoville et al. 2001), which has 
significantly higher spatial resolution than the typical size of an 
H ii region in the luminosity range of interest here. Under these 
circumstances, it is possible to classify a single H n region, due to 
a single coeval OB association, artificially as a set of aggregated 
H ii regions, each of a lower luminosity. 

We thus conclude that blending should not have serious con- 
sequences for the main results presented in this letter, although 
for individual galaxies at larger distances blending may intro- 
duce uncertainties in the LFs. 

3. Calibration 

The continuum subtraction and photometric calibration of the 
Ha images has been described in detail by Knapen et al. (2004). 
We estimated the uncertainties in the calibration to be L = 



0.1 dex by assuming errors of 3% in the distance to the galaxy, 
and of one standard deviation of the adopted background value 
of the calibration star. 

Since the final continuum-subtracted images of most galax- 
ies contain emission from both the Ha line and the [Nn] lines 
at 6548 A and 6583 A, we need to correct for the latter before 
determining the star formation rates. We do this by assuming a 
fixed ratio between the [Nn] 6583 and 6548 lines of three, and a 
value j] for the [Nn] 6583/Ha ratio of 0.25 for galaxies that do 
not contain any or only a few "giant" Hn regions (e.g., early- 
type galaxies) and of 0.16 for the later Hubble types, based on 
descriptions in Osterbrock (1989). We compute the final correc- 
tion factor for each galaxy by additionally determining the rela- 
tive filter transmissions in the Ha and [Nn] lines, for which we 
used the filter transmission characteristics as measured by the 
Isaac Newton Group (ING), which operated the JKT. 

We found that the [Nn] contribution varied due to filters 
used, recessional velocity, and galaxy type, and ranged from 9% 
to 36% with a mean of 17%. NGC 3631 and NGC 4321 were 
both observed with Ha filters of width 15 A which was narrow 
enough not to include any [Nn] emission. 

Thilker et al. (2002) use values of rj which vary from 0. 1 to 
0.4, but only one galaxy is common to both their and our sample. 
This is NGC 5457, for which Thilker et al.'s value is 0.2, as 
opposed to ours of 0.16. Using the former value would lead to a 
difference of 3% in the [Nn] contribution to our image, which is 
not a significant factor in the analysis of the LF. No other such 
comparisons could be made from the literature. 

4. H ii region catalogues 

We used a semi-automatic method to derive the H n region cata- 
logues for all sample galaxies which uses the region programme 
(C. H. Heller, private communication; see Rozas et al. 1999; 
Pleuss et al. 2000). Using a tree algorithm, region catalogues 
each H n region, yielding the position of its centre, its area in 
pixels and its background-subtracted flux. An advantage of re- 
gion over any fully automated programme is that it allows for for 
defining local backgrounds, and for manual editing of the result- 
ing catalogue to remove features like cosmic rays or remnants of 
imperfect continuum subtraction of foreground stars. We chose 
not to catalogue H ii regions in the central 1 kpc region of each 
galaxy for two reasons. Firstly, crowding and blending is espe- 
cially severe there, and secondly, especially in the AGN hosts in 
our sample, the Ha can be severely contaminated by emission re- 
sulting from shocks or AGN emission. The H n region catalogues 
are all available electronically through the Centre de Donnees 
Stellaires (CDS) 1 . The catalogues list for all galaxies a number 
identifying the Hn region (col. 1), the position of the pixel of 
maximum intensity of the Hn region relative to the centre of 
the galaxy, in RA and dec (pixel scale is 0.241 arcsec pixeL 1 ; 
cols. 2 and 3), the logarithm of its integrated luminosity in erg s _1 
(col. 4), the area of the H n region in units of pc 2 (col. 5), and the 
value of the local background (col. 6). 

5. Luminosity functions 

5. 1 . Individual galaxies 

On the basis of our H n region catalogues, we constructed lumi- 
nosity functions (LFs) for all our galaxies, which we present in 

1 Via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 
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Fig. 2. Hn region LFs of NGC 3631 (top panel), whose cata- 
logue contains 801 Hn regions, fitted with a double power law, 
and of NGC 1068 (lower panel), with 166 Hn regions, and fit- 
ted with a single power law. The LFs have been plotted with bin 
widths of 0.2 and 0.5 dex, repectively. 



Fig. 1 (online only). The bin widths for plotting and fitting were 
determined following Scott (1979), and vary from 0.2 to 0.7 dex 
as a function of the number of Hn regions (see Table[Q. We fit 
slopes to the LFs, following the equation N(L)dL = ALr"dL, 
where N{L)dL is the number of H n regions with luminosities 
in the range L to L + dL, and a is the slope, as introduced by 
KEH. Consistent with KEH and most of the more recent work, 
we made an adjustment of -1 to the slopes of the plotted LFs, 
which is due to the fact that we show differential LFs as deter- 
mined with logarithmic binning, whereas the slopes refer to a 
differential LF with linear binning. The slopes, as determined 
with a weighted fit, have been plotted on the LFs in Fig. 1, and 
are listed, along with their formal errors and coefficients of de- 
termination r 2 , in Table 1 . 

KEH pointed out (see also, e.g., Rand 1992; Beckman et al. 
2000) that for those galaxies which have large enough numbers 
of H ii regions with L > 39 dex, the LFs were, in general, best 
fitted with a double power law, with the fit to the higher- L portion 
of the LF steeper-sloped than that of the lower-L portion. This is 
borne out in the present sample, where only eight of our galaxies 
satisfy the requirement of sufficient high-L H n regions to entail 
a statistically significant double power law fit. Fitting the LFs 
of these galaxies with a single power law gives a demonstrably 
inferior fit. Parameters of the single power law fits are listed in 
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Fig. 3. LF made from the combined Hn region catalogues of 53 
of our sample galaxies. The large number of Hn regions allows 
plotting the LF in steps of 0. 1 dex. The best double-slope fit is 
indicated by the lines, where dashes show extrapolations beyond 
the range over which the best fit has been made. 



Table 1 for all those galaxies in the sample with sufficient H n 
regions to do this reliably, whereas the results of the double fits 
are listed in Table 2 (online only). Fig. 2 shows representative 
examples of a galaxy where a double fit is superior {top panel) 
and of one with insufficient high-L H n regions (lower panel). 

5.2. Composite luminosity functions 

In the present Letter, our purpose is to discuss the combined 
LF for the complete sample (rather than to use separately the 
individual LFs shown in Fig. 1). We can produce a combined 
LF because all our data have been uniformly calibrated pho- 
tometrically, and because, separating the sample into three ba- 
sic morphological blocks: early, intermediate and late types, 
we did not find significant differences in the slopes of the low 
or high luminosity portions of the LF. The result, which con- 
tains 17,797 Hn regions, is shown in Fig.|3](bin width 0.1 dex). 
Overlaid on that Figure is the best double power law fit (a U pper = 
-1.36 ± 0.02, r 2 = 0.98; ai ower = -1 .86 + 0.03, r 2 = 1 .00). A sin- 
gle power law can be fitted (a s i ng ie = -1.58+0.01, r 2 = 0.94), but 
has a significantly lower coefficient of determination. The low-L 
limit to the fitted range occurs at a value of L — 37.6 dex, which 
is a completeness limit imposed by a combination of signal-to- 
noise ratio and spatial resolution in the Ho- images. As we will 
see below, this completeness limit occurs well below the break 
luminosity and allows fits to both the low and high-L slopes. 

Not only does the statistical strength of our combined LF of 
almost 18,000 Hn regions allow us to confirm so clearly that 
a break occurs in the LF slope, but the quality of our data is 
good enough for us to establish exactly where that break occurs. 
We made successive fits to the two slopes, varying the break lu- 
minosity between L = 37.7 dex and 39.1 dex. Fig.0]shows the 
parameter which best indicates the quality of each fit, the coef- 
ficient of determination, r 2 , as a function of log L for both the 
lower and upper slopes. Before interpreting the results, we must 
note that r 2 will move closer to unity (indicating better fits) as 
more points are fitted, which will occur for higher break-L for 
the lower-L slope, and vice versa for the upper slope. Even so, 
Fig.|4]clearly shows a maximum for each slope, indicating a dis- 
tinct and restricted range of log L where the best fit quality is 
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Fig. 4. Coefficient of determination, r 2 , for the upper and lower 
slopes of the double-slope fit to the combined LF of Fig. [5] as a 
function of the log L at which the break between the two slopes 
occurs. Better fits are closer to r 2 = 1, so the best fit is achieved 
for a break at log L(erg s _1 ) = 38.6 + 0.1, the average of the best 
values for the low and high-L slopes. 



obtained. The resulting value for the break, or Stromgren, lumi- 
nosity is log L(erg s _1 ) = 38.6 + 0.1. 

Theoretical explanations for the break have been discussed 
previously, notably by Beckman et al. (2000) who suggested that 
it might be caused by the transition from ionisation bounding to 
density bounding as the H n regions become more luminous, and 
who made semi-quantitative predictions of the relations between 
the slopes of the two power laws based on this scenario. Pleuss 
et al. (2000) argued that the effect might be due to a cluster- 
ing property in which H n regions of increasing luminosity have 
an increasing probability of absorbing smaller Hn regions by 
overlap, and that this may be subject to a scaling law related to 
the separation between star formation regions. Any quantitative 
attempt to implement the latter or similar proposals, however, 
leads to an LF in which the upper luminosity slope is shallower, 
not steeper, than that in the lower luminosity range. So the den- 
sity bounding scenario retains its basic plausibility, although it 
must be carefully quantified to take into account the fact that H n 
regions are highly clumpy (see, e.g., Giammanco et al. 2004 for 
a study of how dumpiness affects the propagation of ionising 
radiation in H n regions), rather than being idealised Stromgren 
spheres of uniform density. An alternative explanation is the one 
proposed by Oey & Clarke (1998), in which the break in the LF 
is due to the evolution of the Ha luminosities of the H n regions. 

It is interesting to note that at the high luminosity end of 
the LF, above L = 40dex, the curve falls away quite sharply. 
It is true that the number of Hn regions per bin in this range 
is small, and this is shown clearly in Fig. [3] in the increasing 
amplitude of the error bars as L increases. The LF would take 
this form if we were seeing a true cut-off in the luminosity of H n 
regions due to a physical limit on the masses of high mass stars 
modulated by a statistical effect due to the variation of the total 
stellar mass in young clusters. Before reaching a firm conclusion 
on this, however, it would be useful to be able to increase the 
statistical base in this range by further observations. 



6. Summary and concluding remarks 

We have analysed continuum-subtracted Ha images of 57 
nearby galaxies, and here present the LFs and catalogues of their 
populations of H n regions. Using the combined data set we have 
produced a single LF covering the range in Ha luminosity from 
L = 37.5 dex to just over L = 40dex. The total number of Hn 
regions contributing to this LF is a little under 18,000, which is 
almost an order of magnitude greater than that of any previously 
published LF of this type. The best functional fit of the LF, above 
the completeness limit, is by a double power law, having a break 
at L — 38.6 dex, with the steeper slope to higher luminosity and 
the shallower slope to lower luminosity. The existence of this 
break suggests a change in physical regime going from lower to 
higher Ha luminosity. If confirmed, and by virtue of its occur- 
rence at high H n region luminosities, the break is of potential 
use as a distance indicator for star-forming disk galaxies. 

The number of H n regions in our sample is large enough 
to suggest that the steep fall off in the LF above L = 40 dex 
indicates that we are sampling the range where there is a cut-off 
in the luminosities of the H n regions, implying a drop-off in the 
LF of the underlying massive stellar clusters. This is giving us 
potentially important information about a physical limit to the 
masses of stars in young clusters of varying total stellar mass. 
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Table 2. (Online only) Double fits to LFs for those galaxies 
whose LFs allowed such fits. 
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Fig. 1. Online only LFs of 56 of our sample galaxies, with best 
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Table 1. (Online only) General properties of the sample galaxies, Ha images and resulting Hn region catalogues. Calibration 
constant, in log of L (erg s~') per count in the Ha image (col. 2), lower completeness limit of the LF (log L; col. 3), number of H ii 
regions catalogued (col. 4), image scale in pc per arcsec (col. 5), distance to the galaxy in Mpc (from Paper II; col. 6), and resolution 
(seeing, FWHM) in the Ha image in pc (col. 7). N/A in col. 2 means no calibration available, in col. 3, PQ means poor quality LF 
and NDR means no disk Hn regions. Cols. 8-9 tabulate the upper and lower limits (in logL) of the luminosity function; col. 10 the 
bin width used for the plotting and fitting of the LF, ranging from 0.2 to 0.7 dex as a function of the number of H n regions; col. 1 1 
the value of the slope, a, of a single power law fitted to the data; col. 12 the combined error (in log L), including background effects 
and calibration errors, and col. 13 the coefficient of determination, r 2 , where values closer to unity indicate better fits. 
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